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Abstract. We present high spectral-resolution optical spectra of 49 Herbig Ae/Be stars in a search for the [O i] 6300A line. 
The vast majority of the stars in our sample show narrow (FWHM < 100 km s -1 ) emission lines, centered on the stellar radial 
velocity. In only three sources is the feature much broader (~ 400 km s ), and strongly blueshifted (-200 km s ) compared to 
the stellar radial velocity. Some stars in our sample show double-peaked lines profiles, with peak-to-peak separations of ~ 10 
km s . The presence and strength of the [O i] line emission appears to be correlated with the far-infrared energy distribution of 
each source: stars with a strong excess at 60 yum have in general stronger [O i] emission than stars with weaker 60 jim excesses. 
We interpret these narrow [Oi] 6300A line profiles as arising in the surface layers of the protoplanetary disks surrounding 
Herbig Ae/Be stars. A simple model for [O i] 6300A line emission due to the photodissociation of OH molecules shows that 
our results are in quantitative agreement with that expected from the emission of a flared disk if the fractional OH abundance is 
~5 x 10- 7 . 
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1. Introduction 

Herbig Ae/Be (HAEBE) stars are intermediate-mass pre-main- 
sequence objects. The spectral energy distribution (SED) of 
these sources is characterized by the presence of an infrared 
excess, due to thermal emission of circumstellar dust. For 
HAEBE stars of spectral type F, A and late-B, the evi- 
dence for a disk-like geometry of this circumstellar dust (e.g . 
Mannings & Sargentlll997tlTesti et alJkOoH IPietu et alJ fcoOS ; 
Fuente et al. 2003; Natta et al. 2004) is generally accepted. The 
spatial distribution of the circumstellar matter a round early-B- 
type stars is less clear (e.g. iNatta et aljEo OO). Early-B stars 
have dissipation time scales for the circumstellar spherical 
envelope of the order of their pre-main-sequence life time. 
Nevertheless observation s seem to indicate that at least some of 
these sources have disks dVink et a l. 2002). Both disk-like and 
spherical geometries might be present around these sources. 



* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (program numbers 54.D-0363 and 68. C- 
0348). 

** Visiting Astronomer, Kitt Peak National Observatory, National 
Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under coop- 
erative agreement with the National Science Foundation. 



iMeeus et al] feOOll henceforth M01) have classified the 14 
isolated HAEBE stars in their sample based on the shape of 
the mid-IR (20-100 pm) SED. Group I contains the sources 
which have a rising mid-IR flux excess (double-peaked SED). 
Group II sources display more modest IR excesses. M01 sug- 
gest phenomenologically that the difference in SED shape re- 
flects a different disk geometry: group I sources have flared 
disks, while group II objects have geometrically flat disks. 

iDullemondl J2002I) IDullemond et al J J2002I) and 
iDullemond & Dominikl d2004l) have modeled circumstel- 
lar disks with a self-consistent model based on 2D-radiative 
transfer coupled to the equations of vertical hydrostatics. The 
model is composed of a disk with an inner hole (~ 0.5 AU), 
a puffed-up inner rim and an outer part. The inner rim, which 
is located at the dust sublimation temperature, is puffed-up 
due to the direct head-on stellar radiation and the consequent 
increase of the local gas temperature. The models show that 
th e outer part o f the disk can either be flared (as in the models 
of IChiang & Goldreichl 19971) . or lie completely in the shadow 
of the puffed-up inner rim. The SED of a flared disk displays a 
strong mid-IR excess, comparable to that of the M01 group I 
sources. Self-shadowed disk models have SEDs with a steeper 
decline towar ds longer wa velengths, like the SEDs of group II 
sources. The IDullemondl d2002l) models hence connect the 
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M01 classification to the disk geometry: group I sources have 
fla red disks, group II objects self-s hadowed disks. 

lAcke & van den Anckeil J2004 henceforth AV04) have 
investigated all a vailable near-infra red Infrared Space 
Observatory (ISO, iKessler et al] 1 1996b spectra of HAEBE 
stars. In their paper, the emission bands at 3.3, 6.2, "7.7", 
8.6 and 11.3 fim were studied qualitatively. These emission 
features are attributed to po l ycyclic aromatic hydroca rbons 
(PAHs, iLeger & Pugetl Il984 lAllamandola et alJ Il989l) . The 
latter, carbonaceous molecules, are thought to be excited by 
ultraviolet (UV) photons, and radiate in the near-infrared 
in bands linked to the vibrational modes of the CC and CH 
bonds. AV04 have found a correlation between the strength of 
the PAH features and the shape of the SED: group I sources 
display significant PAH emission in their spectra. Group II 
members on the other hand show much weaker PAH bands, 
indicating the lack of excited PAH molecules in these systems. 
Following the remark of M01, AV04 suggest this is due to 
the geometry of the circumstellar disk. In self-shadowed 
geometries, the stellar UV flux needed to excite the PAH 
molecules cannot reach the disk surface, while in flared disks 
this surface is directly exposed to the radiation field of the 
central star. Theoretical modeli ng of the P AH emission in 
flared circumstellar disks jHabartetalj r2n04) is in agreement 
with the observations. 

Forbidden-line emission from different element s (C i, N 
n, O i, O ii, S ii, Ca n, Cr n, Fe n, Ni n; see e.g., iHamannl 
1994) in the optical part of the spectrum has been observed 
in a considerable amount of HAEBEs. In the present paper 
we focus on the [Oi] features at 6300 and 6363 A . The cir- 
cumstellar region from which these lines emanate has been 
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blueshifted high-velocity (a few 100 km s _1 ) wing observed 
in a few sources is formed in an outflow whose redshifted 
part is obscured by the circumstellar disk. The origin of 
low- velocity symmetric [ Oil emission profile is less clear. 
iKwan & T ademarui Jl988l) have provided a qualitative explana- 
tion for the observed [Oi] profiles in T Tauri stars. Their model 
consists of a low- and a high-velocity component, the first due 
to a slow disk wind, the latter emanating from a collimated jet. 
Previously published suggestions for the forbidden-line emit- 
ting region in HAEBE stars are based on this m odel and include 
a low -velocity disk wind (iHirfh et alJll994 ICorcoran & Ravi 
1 19971) and a spherically symmetric stellar wind {without an ob- 
scuring circumstellar disk, IBohm & Cataialll994 . In the latter 
explanation, the authors suggest that the [0 1] line is formed in 
the outermost parts of the stellar wind. 

With the present paper, we intend to contribute to this dis- 
cussion, and try to explain the observed low-velocity compo- 
nent in a broader framework. In |2] the composition of the 
sample is presented. First, we investigate the observational data 
and perform a quantitative analysis. We search for correla- 
tions between parameters describing the forbidden-line emis- 
sion and the SED. We have also checked the possible connec- 
tion between the presence and strength of PAH features and the 



forbidden-line emission (§g3j|4|l. Second, we propose a simple 
model for the forbidden-line emission region, and compare the 
model results to the observations (|5j. 

2. The data set 

2.1. The sample 

We have observed 49 Herbig Ae/Be stars in the wavelength 
region around the forbidden oxygen line at 6300A . Since we 
intend to compare the results of the ISO study of HAEBE stars 
by AV04 to the present [Oi] -emission-line analysis, we have 
compiled our [0 1] sample in order to have as large as possible 
an overlap with the ir sample. Additional HAEBE stars from the 
catalogues oflThe et alJJl994 Table 1 and 2 in their article) and 
Mal fait et alJ ( 119981) were observed when possible to enlarge 
the sample. 

Optical spectra of the sample stars were obtained with four 
different instruments on five different telescopes: the Coude 
Echelle Spectrometer (CES) on the Coude Auxiliary Telescope 
(CAT) and on the ESO 3.6m telescope, the Fiber-fed Extended 
Range Optical Spectrograph (FEROS) on the ESO 1.5m tele- 
scope, the CCD Echelle Spectrograph on the Mayall 4m tele- 
scope at Kitt Peak National Observatory (KPNO) and the 
Utrecht Echelle Spectrograph (UES) on the William Herschel 
Telescope (WHT). The first three telescopes are located at the 
ESO site La Silla (Chile), KPNO is in Arizona (USA) and the 
WHT is situated on La Palma (Canary Islands, Spain). The 
FEROS and WHT measurements were made by Gwendolyn 
Meeus. The spectra were obtained at different periods during 
the past 10 years. In Tabled we have summarized the sample 
stars and spectra included in this analysis. For 10 sources, two 
or more spectra were obtained. 

In order to extract the spectra from the raw data, a stan- 
dard echelle-data reduction was applied. This includes back- 
ground subtraction, cosmic-hit removal, flatfielding and wave- 
length calibration. The spectra were normalized to unity by fit- 
ting a spline function through continuum points and dividing 
by it. The reduction and normalization was performed using 
the ESO program MIDAS. 

The spectral resolution A/AA at 6300A is instrument- 
dependent. The resolution for the KPNO data is ~30,000, for 
the FEROS and WHT spectra -45,000, and -65,000 for the 
CAT data. The highest spectral resolution is reached in the ESO 
3.6m spectra with A/AA * 125,000. The FEROS and WHT 
data cover a large fraction of the optical wavelength range 
(3700-8860A and 5220-91 10A respectively), which allows us 
to accurately determine the radial velocity of the central star 
(vmd) based on many photospheric absorption lines through- 
out the spectrum. The KPNO data cover the wavelength range 
between 5630 and 6640A. We have used the Fe n absorption 
line at 5780. 128 A to estimate the stellar radial velocity in these 
spectra. The CES spectra (CAT and ESO 3.6m) are limited 
to one order around the [Oi] line at 6300A. No stellar radial 
velocities could be determined from these spectra. For a few 
sources, Ca 11 K 3934A spectra are at our disposal. When no 
photospheric-line spectra were available to us, we estimated 
the radial velocity based on these circumstellar lines. An esti- 



Acke, van den Ancker & Dullemond: [O i] emission in HAEBE disks. 



3 



mate of v rac i was retrieved from the literature in case we could 
not determine it from any of our spectra. When no literature 
estimate exists, we have computed an average v ra d of stars in 
the area on the sky (radius 2') around the source. The radial 
velocities are included in Table|4](see later). 

The spectra were velocity-rebinned and, after applying a 
heliocentric correction, centered around the stellar radial ve- 
locity. In this way we can compare the measured velocities, 
independently of the intrinsic radial velocity of the entire sys- 
tem. 

In this analysis, we focus on the [Oi] emission line at 
6300A. For the sample stars with FEROS, KPNO and WHT 
spectra, the extended spectral range allows us to include the 
[Oi] line at 6363A as well. Furthermore, these spectra contain 
the Ha line at 6550A. 

2.2. Measurement of the [Oi] emission lines 

Before analysing the [O i] 6300A line, we have removed the 
telluric absorption lines from the spectrum by cutting out these 
wavelength regions and replacing them by a spline approxima- 
tion based on all the other data points in the spectrum. The air- 
glow feature ([Oi] emission from the Earth's atmosphere) was 
suppressed in the same way. 

In ten of our sample sources, the 6300A region is rich in 
photospheric absorption lines. This hampers the detection of 
possible superimposed [O i] emission. Of the other 39 stars, 29 
objects display a pure-emission shape (i.e. no underlying ab- 
sorption features). The ten remaining sources have no detected 
[Oi] line. 

We have measured the equivalent width (EW) in A and cen- 
troid position 1 in km s _! . Following tradition, the EW is nega- 
tive when the feature is in emission. The EW in absorption-rich 
6300A spectra is determined in a fixed interval (±50 km s _1 ) 
around the stellar radial velocity. We remind the reader that 
the centroid position is given with respect to the reference 
frame of the central star: a positive (negative) centroid posi- 
tion indicates that the feature is redshifted (blueshifted) com- 
pared to the stellar radial velocity. For the detected emission 
lines, we have also determined the full width at half maxi- 
mum (FWHM) by fitting a Gaussian function to the feature. 
This estimate of the FWHM is corrupted by the instrumen- 
tal profile of the spectrograph. Since the telluric lines in the 
spectra are intrinsicly very narrow, the FWHM of these ab- 
sorption features is a good estimate for the instrumental width. 
The measured FWHM of the [O i] lines is corrected for this in- 
strumental broadening. All detected features in this sample are 
spectrally resolved (i.e. have a corrected FWHM several times 
larger than the instrumental width). Finally, we have defined 
an asymmetry parameter J[ based on the centroid position and 
the extreme blue (S) and red (K) ends of the emission profile: 
£F{ = (centroid -3) ICR- centroid). When J[ is larger (smaller) 
than 1, the centroid position lies closer to the red (blue) end of 
the emission line, hence indicating a stronger red (blue) wing 

1 centroid position = j'/"^^! with la the continuum-subtracted in- 
tensity of the profile and v the velocity parameter. 



compared to the center C of the line 2 . Note that the line profile 
can be completely blueshifted compared to the stellar radial ve- 
locity, while having The only [O i] line parameter that is 
sensitive to (errors on) the stellar radial velocity is the centroid 
position. 

Errors on the [O i] parameter values have been estimated 
based on the noise in the spectrum. In the figures that are pre- 
sented in this paper, we have chosen to show a representative 
error bar instead of plotting the individual error bars on each 
measurement. This was done for the sake of clarity. The plot- 
ted error bar, which indicates the errors on a mean entry in the 
figure, is refered to in the captions of the figures as the typical 
error bar. 

When no [O i] feature was detected, we have computed an 
upper limit for the EW from the noise on the data (cr). We 
assume that the strongest feature that remains undetected is a 
rectangular emission line with a height of 5 times the noise on 
the data and a width of 1 A (= 47 km s ). The equivalent width 
of this hypothetical feature thus is EW,„ flA = 5cr A. The latter 
value was used as the upper limit for the undetected line. 

For the 10 sources for which we have more than one opti- 
cal spectrum at our disposal, we have measured the [O i] 6300A 
feature in all spectra separately. Afterwards the parameter val- 
ues were compared. No significant differences were noted, ex- 
cept for Z CMa (for which the |EW| increased with ~40% in 
6 years) and MWC 297 (decrease of |EW| with -25% in 3 
months). This suggests that for most stars in our sample the 
[O i] emission is constant in time, although variations -either 
in the continuum flux or in the line emission itself- do occur 
in some objects. The final data consist of one measurement for 
each source, which is a weighted average when two or more 
spectra are available. The centroid positions of the features in 
the different spectra of the same source agree well within the 
error bar, except again for MWC 297, where the difference 
in centroid position adds up to 20 km s _1 . The forbidden-line 
emission region of this source apparently rotates as a whole 
around the central star. 

Thanks to the large wavelength coverage of the FEROS, 
KPNO and WHT spectra, also the wavelength region around 
the forbidden [O i] 6363 A emission line was observed for 30 
objects. When present in these spectra, we have determined the 
EW and FWHM of the [O i] 6363 A emission line. Otherwise, 
an upper limit for the EW was determined. The EW of the Ha 
line was measured as well. When no Ha spectrum was avail- 
able to us, we have included literature values. 

2.3. The SED parameters 

To characterize the spectral energy distribution (SED) of the 
central star, several quantities were determined, based on UV- 
t o-millimeter p hotometry from the literature. For each source, 
a Kurucz (1991) model with effective temperature T e ff and sur- 
face gravity \ogg corresponding to the spectral type of the 
source was fitted to the de-reddened photometry. From this 
model, parameters like total luminosity L and UV luminosity 
L uv (2-13.6 eV) can be determined. By subtracting the Kurucz 

2 C = (« + S)/2 
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Table 1. The sample of HAEBEs used in this analysis, based on iThe et alJ i 19941) and iMalfait et alJ dl998h . For each object, 
the available optical spectra are indicated. Furthermore, the SIMBAD values for right ascension (RA) and declination (Dec) are 
given, as well as the date, start time and integration time T of the observation. 



Sample stars 


Object 


Spectrum 


RA (2000) 


Dec (2000) 


Date 


Start 


T 






h m s 


o r rr 


dd/mm/yy 


h m 


[m] 


V376 Cas 


KPNO 


00 11 26.1 


+58 50 04 


25/06/02 


10:22 


20 


VXCas 


KPNO 


00 31 30.7 


+61 58 51 


24/06/02 


10:06 


25 


AB Aur 


FEROS 


04 55 45.8 


+30 33 04 


17/01/99 


22:37 


21 


HD 31648 


WHT 


04 58 46.3 


+29 50 37 


24/12/96 


00:00 


20 


HD 34282 


ESO 3.6m 


05 16 00.5 


-09 48 35 


01/04/02 


00:49 


50 


HD 34700 


ESO 3.6m 


05 19 41.4 


+05 38 43 


01/04/02 


00:05 


40 


HD 35929 


CAT 


05 27 42.8 


-08 19 38 


13/01/94 


05:27 


30 


HD 36112 


WHT 


05 30 27.5 


+25 19 57 


24/12/96 


02:18 


25 


HD 244604 


WHT 


05 31 57.3 


+ 11 17 41 


24/12/96 


02:48 


20 


HD 245185 


WHT 


05 35 09.6 


+ 10 01 52 


24/12/96 


03:42 


12 


V586 Ori 


CAT 


05 36 59.3 


-06 09 16 


13/01/94 


03:10 


60 


BFOri 


CAT 


05 37 13.3 


-06 35 01 


13/01/94 


04:22 


60 


ZCMa 


ESO 3.6m 


07 03 43.2 


-11 33 06 


01/04/02 


01:57 


30 




WHT 






24/12/96 


05:36 


20 


HD 95881 


ESO 3.6m 


11 01 57.6 


-71 30 48 


01/04/02 


23:46 


17 


HD 97048 


ESO 3.6m 


11 08 03.3 


-77 39 18 


01/04/02 


22:56 


20 




FEROS 






28/01/99 


04:28 


55 


HD 98922 


CAT 


11 22 31.7 


-53 22 12 


13/01/94 


08:40 


15 


HD 100453 


ESO 3.6m 


11 33 05.6 


-54 19 29 


01/04/02 


01:44 


10 




FEROS 






26/01/99 


04:55 


20 


HD 100546 


CAT 


11 33 25.4 


-70 11 41 


13/01/94 


07:55 


10 




ESO 3.6m 






01/04/02 


02:32 


10 




ESO 3.6m 






31/03/02 


23:27 


10 




FEROS 






28/01/99 


05:27 


30 


HD 101412 


CAT 


11 39 44.5 


-60 10 28 


13/01/94 


09:06 


20 


HD 104237 


ESO 3.6m 


12 00 05.1 


-78 11 35 


31/03/02 


23:35 


5 




FEROS 






27/01/99 


03:38 


20 


HD 135344 


ESO 3.6m 


15 15 48.4 


-37 09 16 


01/04/02 


02:40 


23 




FEROS 






27/01/99 


04:02 


45 


HD 139614 


ESO 3.6m 


15 40 46.4 


-42 29 54 


01/04/02 


03:06 


17 




FEROS 






27/01/99 


04:49 


30 


HD 141569 


ESO 3.6m 


15 49 57.8 


-03 55 16 


01/04/02 


04:55 


6 




KPNO 






23/06/02 


03:57 


10 


HD 142666 


ESO 3.6m 


15 56 40.0 


-22 01 40 


01/04/02 


03:52 


29 


HD 142527 


ESO 3.6m 


15 56 41.9 


-42 19 23 


01/04/02 


03:30 


20 


HD 144432 


ESO 3.6m 


16 06 58.0 


-27 43 10 


01/04/02 


04:22 


23 


HR 5999 


ESO 3.6m 


16 08 34.3 


-39 06 18 


01/04/02 


04:46 


6 


HD 150193 


ESO 3.6m 


16 40 17.9 


-23 53 45 


01/04/02 


05:04 


30 


AK Sco 


ESO 3.6m 


16 54 44.9 


-36 53 19 


01/04/02 


05:34 


29 


CD-42° 11721 


ESO 3.6m 


16 59 06.8 


-42 42 08 


01/04/02 


06:05 


60 


HD 163296 


ESO 3.6m 


17 56 21.3 


-21 57 22 


01/04/02 


07:08 


5 


HD 169142 


ESO 3.6m 


18 24 29.8 


-29 46 49 


01/04/02 


07:15 


17 


MWC 297 


ESO 3.6m 


18 27 39.6 


-03 49 52 


01/04/02 


08:43 


60 




KPNO 






25/06/02 


05:40 


20 


VV Ser 


KPNO 


18 28 49.0 


+00 08 39 


23/06/02 


04:49 


30 


RCrA 


CAT 


19 01 53.7 


-36 57 08 


19/08/94 


12:43 


45 


TCra 


ESO 3.6m 


19 01 58.8 


-36 57 49 


01/04/02 


07:39 


40 


HD 179218 


ESO 3.6m 


19 11 11.6 


+ 15 47 16 


01/04/02 


09:47 


8 




KPNO 






24/06/02 


04:24 


10 


WW Vul 


KPNO 


19 25 58.8 


+21 12 31 


23/06/02 


06:31 


30 


HD 190073 


KPNO 


20 03 02.5 


+05 44 17 


24/06/02 


05:01 


10 


BD+40°4124 


KPNO 


20 20 28.3 


+41 21 52 


23/06/02 


08:11 


20 


LkHff 224 


KPNO 


20 20 29.3 


+41 21 28 


24/06/02 


05:40 


30 


LkHa 225 


KPNO 


20 20 30.5 


+41 21 27 


23/06/02 


09:21 


30 
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Sample stars (continued) 








Object 


Spectrum 


RA (2000) 


Dec (2000) 


Date 


Start 


T 






h m s 




dd/mmfyy 


h m 


[m] 


PVCep 


KPNO 


20 45 54.3 


+67 57 39 


24/06/02 


07:18 


30 


HD 200775 


KPNO 


21 01 36.9 


+68 09 48 


23/06/02 


09:56 


15 


V645 Cyg 


KPNO 


21 39 58.0 


+50 14 24 


25/06/02 


07:12 


30 


BD+46°3471 


KPNO 


21 52 34.1 


+47 13 45 


25/06/02 


08:52 


20 


SV Cep 


KPNO 


22 21 33.3 


+73 40 18 


23/06/02 


10:19 


20 


ILCep 


KPNO 


22 53 15.6 


+62 08 45 


25/06/02 


09:58 


20 


MWC 1080 


KPNO 


23 17 25.6 


+60 50 43 


24/06/02 


08:59 


30 



model from the infrared photometry, the flux excesses at 2.2 
(i.e. K band), 60, 850 and 1300 micron can be determined. 
Other parameters used in this analysis are the the observed 
bolometric luminosity Lboi (not corrected for extinction) and 
the distance d to the source. 

To convert the measured EW of the [O i] lines to [O i] lumi- 
nosities, we have computed the theoretical 6300A continuum 
flux from the Kurucz model. The relationship between the [O i] 
luminosity L([Oi]) and the EW is 

L([Ol]) = 4;rd 2 2.5 x 10 2 |EW| x ^"(6300A) (1) 

where L([Oi]) is in [L Q ], d in [pc], EW in [A] and the contin- 
uum flux T in [W m~ 2 //m -1 ]. In a similar way, Ha luminosities 
were computed from the Ha EWs. 

The sample sources were classified into the M01 groups, 
based on the shape of their IR excess. In accordance with 
Ivan Boekel et alJ J2003I) . we characterize the SED using the 
following quantities: the ratio of L{NIR) (the integrated lumi- 
nosity derived from broad-band J, H, K, L and M photome- 
try) and L(IR) (the corresponding quantity derived from IRAS 
12, 25 and 60 micron points), and the non-color-corrected 
IRAS [12] - [60] color. These quantities naturally separate 
sources with a strong mid-IR excess (group I) from more 
modest mid-IR emitters (group II). In Fig. ^ the classifica- 
tion is visualized in a diagram. The dashed line represents 
L(NIR)/L(IR) = ([12] - [60]) + 0.9, which empirically pro- 
vides the best separation between bot h groups. This classi fica- 
tion method has also been used in Dul lemond et alJ (|2003) and 
AV04 and has been described more thoroughly in lAcke e t al. 
fe004 . The quantities L(NIR)/L(IR) and [12] - [60] are listed 
in Table|2] 

Five sources in this sample display the amorphous 10 mi- 
cron silicate feature in absorption (AV04). Considering the 
presence of an outflow and CO bandhead emission, strong veil- 
ing and a UV excess, these objects are thought to have disks 
whose luminosity is dominated by viscous dissipation of en- 
ergy due to accretion, and are deeply embedded sources. This 
idea is supported by their very red IR spectral energy distri- 
bution. They are likely fundamentally different from the other 
sample stars. We therefore classify them in a different group: 
group III. 

BD+40°4124, R CrA and LkHa 224 have not been classi- 
fied based on their location in the diagram in Fig.^ Confusion 
with background sources in the IRAS photometry made it im- 
possible to derive the desired quantities L(NIR)/L(IR) and 
[12] - [60]. BD+40°4124 has been classified as a group I 
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Fig. 1. Classification diagram, based on Ivan Boekel et alJ 
( 2003). Objects in the upper left corner are classified as group II 
sources (open diamonds), sources on the lower right side are ei- 
ther group I (filled diamonds) or group III (squares) members. 
The dashed line represents the empirical separation between 
group I and II: L(NIR)/L(IR) = ([12] - [60]) + 0.9 (AV04). 
The quantities L(NIR)/L(IR) and IRAS [12] - [60] are listed 
in TablefJ] 

source, because its SED resembles that of HD 200775. R CrA 
and LkHg 224 a re bo th UX Orionis stars. According to 
Dull emond et alJ (|2003), these sources are group II sources. 
Hence we have classified them as members of the latter group. 

In Figs.[2l[3]and|4] the velocity-rebinned spectra of the [O i] 
6300A emission lines are presented. The first plot contains the 
detected emission profiles. Fig. [3] shows the spectra in which 
the possible detection of [O i] emission is confused by underly- 
ing photospheric absorption lines. Sample sources with unde- 
tected 6300A emission are displayed in Fig.0] In each plot, the 
targets are listed according to group. 

3. Analysis of the emission lines 

3.1. The [Oi] emission lines 

The forbidden oxygen 6300A emission line in Herbig Ae/Be 
stars has been studied in previous articles ( Finkenzeller 1985; 



Bohm & Catala 


1 19941 


Hamann 1994; 


Bohm & Hirth 


1997; 


Corcoran & Rav 


|l997l 


1998; Hernandez et al. 2004). We have 



compared the equivalent widths measured in our spectra to 
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Fig. 2. The detected [Oi] 6300 A emission line spectra. Along the x-axis, the velocity in km s _1 is given, on the y-axis the 
normalized intensity. The horizontal dotted line indicates the continuum level, the vertical line the stellar radial velocity (y ra d = 0). 
The source's name and group (in parentheses) are indicated in the upper left corner of each frame. When more than one spectrum 
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Fig. 3. The 6300A spectra of the sample stars for which the detection of [O i] emission is confused by photospheric absorption 
lines. Legend see Fig- EJ 
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Fig. 4. The 6300A spectra of the sample stars with undetected [Oi] emission. Legend see Fig. El 



published literature values in Fig. [5] The error on our EWs is 
computed from the noise on the data, multiplied with the av- 
erage FWHM of the feature. Few previous authors indicate the 
errors on their measurements, hence we cannot set an error bar 
on the literature numbers. Nevertheless we are convinced that 
our data set, consisting of high-resolution, high-signal-to-noise 
spectra (S/N~30-375), are of better quality than the data stud- 
ied in previous papers. We believe that our error bars are a reli- 
able estimate for the scatter on our determinations, and a lower 
limit to the literature errors. The observed scatter in Fig. [5] is 
quite large. This might b e due to intrinsicly va riable [O i] emis- 
sion (e.g. Z CMa. lvan den Ancker et all2004l) or variable con- 
tinuum emission (e.g. in UX Orionis stars). A possible extrin- 
sic explanation is the uncorrected airglow emission in previous 
articles. Since we have removed the telluric lines and airglow 
wavelength regions from our data, our measurements are less 
affected by these features. Low spectral resolution does not al- 



low for removal of the latter. Most of the literature EW val- 
ues in Fig.|5]-especially those of the weak lines- appear to be 
larger than our values, suggesting that the airglow effect indeed 
corrupts these measurements. 

In Section 12.21 we have categorized the profiles based on 
the presence of absorption features underlying the [O i] emis- 
sion line at 6300A. The latter features are photosperic absorp- 
tion lines, which have no direct relation with the [O i] emission. 
This idea is supported by the fact that all 10 sources which dis- 
play absorption features in the 6300A region have late spectral 
types (i.e. low effective temperatures; 6250K < T e fi < 8400K). 
Such objects display photosperic absorption lines in this wave- 
length region, contrary to early-type stars which have flat con- 
tinua around 6300A. 

In Fig. |6] the measured EWs of the line profiles is plotted 
versus the effective temperature. In the left part of the figure, 
the targets with detected [O i] emission are plotted. These ap- 
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Table 2. The non-color-corrected IRAS [12] - [60] color and 
the near-to-mid IR luminosity ratio L(NIR)/L(IR) of the sam- 
ple sources. 



Object 


IRAS [12] - [60] 


L(N1R)IL(1R) 


Group I 






V376 Cas 


1.42 


0.22 


AB Aur 


1.50 


1.73 


HD 34282 


3.04 


2.65 


HD 34700 


3.45 


1.83 


HD 36112 


1.75 


2.08 


HD 245185 


0.24 


0.46 


HD 97048 


1.81 


0.57 


HD 100453 


1.84 


1.06 


HD 100546 


1.00 


0.27 


HD 135344 


2.94 


3.47 


HD 139614 


1.65 


0.88 


HD 142527 


2.51 


1.81 


CD-42° 11721 


3.05 


0.25 


HD 169142 


2.50 


1.00 


TCra 


2.28 


0.42 


HD 179218 


0.27 


0.59 


HD 200775 


2.83 


0.88 


MWC 1080 


1.98 


1.79 


Group II 






VXCas 


-0.12 


0.85 


HD 31648 


0.09 


2.51 


HD 35929 


-1.02 


21.56 


HD 244604 


-0.67 


2.73 


V586 Ori 


-0.29 


2.06 


BFOri 


0.47 


3.23 


HD 95881 


-1.97 


2.69 


HD 98922 


-2.03 


2.42 


HD 101412 


-0.70 


1.80 


HD 104237 


-0.49 


3.55 


HD 141569 


2.50 


8.99 


HD 142666 


-0.16 


1.73 


HD 144432 


-0.29 


2.28 


HR 5999 


-0.85 


5.40 


HD 150193 


-0.96 


1.63 


AKSco 


0.85 


2.63 


HD 163296 


0.41 


2.42 


VV Ser 


-0.39 


2.80 


WW Vul 


-0.10 


2.38 


HD 190073 


-1.44 


3.37 


BD+46°3471 


0.27 


5.39 


SV Cep 


-0.39 


0.68 


IL Cep 


<0.71 


>3.16 


Group III 






ZCMa 


0.87 


0.95 


MWC 297 


1.90 


1.05 


LkHa 225 


1.66 


0.03 


PV Cep 


1.48 


0.37 


V645 Cyg 


1.27 


0.33 



pear to be the early-type sources. For stars in which the possible 
[O i] 6300A emission is confused by photospheric lines (right 
part of the plot), the EWs increase with decreasing . This is 
probably the combined effect of the increasing strength of the 
photospheric absorption lines on one hand, and the decreas- 
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Fig. 5. Comparison between published values for the equiva- 
lent width of the [Oi] 6300A line and our determinations. The 
different plotting symbols refer to data obtained by different 
authors. The typical error bar is plotted in the upper left corner. 
The error bar on the literature values is expected to be larger 
than the indicated error, which is equal to our estimate for the 
error bar. The dashed line represents A EW = 0. 



ing UV luminosity in late-type stars on the other (see later). 
Because of the difficulty to distinguish between the possible 
[O i] emission line and the underlying photospheric absorption 
spectrum, no accurate measurements of the FWHM and cen- 
troid position could be made for these objects. Therefore, these 
determinations have been left out of the further analysis. 

The other sample sources display no absorption lines in the 
[Oi] wavelength range. 10 objects do not display any signifi- 
cant features and will be included in the figures as upper limits. 
The majority of our sample sources (29/49) have a clear [Oi] 
emission profile in their spectra. In Figs.[7]to[K)] the parameters 
that describe the emission [O i] profile are plotted against each 
other. The mean parameters describing the emission line pro- 
file are (FWHM) = 47 km s _1 , (centroid position) = -4 km s _1 
and (&{) = 1 .04. The median values of these quantities in this 
sample are 47 km s _1 , -7 km s'and 0.95 respectively. It is 
striking that, when excluding the three outliers (PV Cep, V645 
Cyg and Z CMa; M in the plots) which will be discussed later, 
most of the sample sources with pure-emission [O i] profiles lie 
close to these average values. The objects seem to have quite 
uniform [Oi] emission line profiles. The EW of the feature is 
not correlated to its width (FWHM), indicating that the broad- 
ening mechanism of the line does not depend on the amount 
of [O i] emission. Nevertheless, a difference in FWHM is ob- 
served between group I (mean: 36 km s , median: 34 km s _I ) 
and group II (mean: 67 km s _1 ; median: 55 km s _1 ). The cen- 
troid positi on is also independent of the EW (Fig|8}. 

iBohm & Caiaial ( 1 1994 note that the [O i] 6300A profiles in 
their sample a re not blue shifted with respect to the stellar ra- 
dial velocity. Icorcoran & Ravi i 19971) on the other hand report 
that most of the HAEBEs in their sample show low-velocity 
blueshifted centroid positions. We look into this more care- 
fully. Fig. |9]displays the distribution of the centroid positions 



Acke, van den Ancker & Dullemond: [O i] emission in HAEBE disks. 



< 



□ 



o 

□ 

>♦ ♦ 



o ♦ 
♦ 

o 



♦ ♦ 



iog(T eff [K]) 

Fig. 6. The EWs of the [O i] emission lines versus the effective temperature of the central star. Filled diamonds refer to group I 
sources, open diamonds to group II. Squares are group III sources, with the x-marks indicating the three sources with an extended 
highly-blueshifted wing. On the y-axis, |EW| is plotted, on the x-axis the logarithm of the effective temperature. The typical error 
bar is indicated in the upper right corner. The right panel contains the 10 sample stars for which the 6300 A region is rich in 
photospheric lines. The EWs are positive due to dominance of the absorption spectrum. The y-axis is flipped to make clear that 
increasing EW designates an increasing contribution of the absorption features. The other 39 sample sources (of which 29 with 
detected [Oi] 6300A emission) are plotted in the left panel. These objects have earlier spectral types. 
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Fig. 7. The EW of the detected [Oi] emission lines versus 
the FWHM. Filled diamonds represent group I, open diamonds 
group II and squares group III sources. The x-marked squares 
represent PV Cep, V645 Cyg and Z CMa. The dash-dotted line 
indicates the mean FWHM (47 km s _1 ). The dotted and dashed 
lines indicate the mean values of group I (34 km s -1 ) and 
group II (67 km s _1 ) respectively. The typical error bar is rep- 
resented on the right side of the plot. 



of the pure-emission [O i] profiles in a histogram. The group III 
sources have been excluded from the diagram. The width of the 
bins is 10 km s _1 . The estimated error on the measured centroid 
positions, including the error on the stellar-radial-velocity de- 
termination, is cr — 15 km s _1 . We assume that the errors on the 
centroid positions are distributed following a standard-normal 
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Fig. 8. The EW of the detected [O i] emission lines versus the 
centroid position of the feature. The plotting symbols are as in 
Fig- El The vast majority of the features is centered around the 
stellar radial velocity within the error bars. The dashed line rep- 
resents the average centroid position (-4 km s _1 ). The typical 
error bar is represented in the lower left corner. 



distribution D exp . We check whether the observed distribution 
Z) bs is representative for the expected normal distribution D exp . 
The possibility that the centroid position of a feature is > cr un- 
der the hypothesis of the D exp distribution is 16%. Five sources 
display redshifted centroid positions larger than 15 km s _1 . 
This is not a statistically significant deviation from D exp for a 
sample of 24 stars. Nevertheless, on the blue side, nine sources 
display centroid positions < -15 km s _1 . The possibility for 
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Fig. 9. The histogram of the centroid positions of the 12 
group I and 12 group II sample sources with a detected [Oi] 
6300 A feature. The bin size is 10 km s . The dashed line rep- 
resents distribution D exp , the full line is the Gaussian fit to the 
observed distribution. The cross-hatched part of the bars indi- 
cate group I sources, the open parts group II members. 



this to happen under the current hypothesis is less than 1%. 
We conclude that, although the majority of stars in our sam- 
ple have centroid positions that are compatible with the stel- 
lar radial velocity, there appear to be some objects that indeed 
show low- velocity blue shifted centroid positions as reported by 
ICorcoran & Ravi i 19971) . The Gaussian fit (full line in Fig. [9} to 
the observed histogram is centered around -7 km s' 1 and has 
cr - 21 km s" 1 . Notice that for the five sources with detected 
[Oi] 6300A emission for which we have spectra covering a 
large part of the optical wavelength region (i.e. AB Aur, HD 
245185, HD 97048, HD 100546 and Z CMa), the error on the 
centroid determination is much smaller than 15 km/s. Of these 
objects, only group III member Z CMa shows evidence for a 
blueshifted feature. In particular the small error (~2 km s _1 ) 
for HD 100546 proves that at least some sources are definitely 
not blueshifted. 

Even though some sources display blueshifted emission, 
there does not seem to be a strong deviation from asymme- 
try. In Fig.E3the FWHM of the [Oi] feature is plotted versus 
the asymmetry parameter M. No correlation between the two 
parameters is noted. 

Three group III sources (PV Cep, V645 Cyg and Z CMa) 
have radically different [O i] emission profiles from the other 
sources in our sample. The shape of these features is double- 
peaked, with one peak close to the radial velocity of the central 
star (v = km s _1 ) and the other at high blueshift (v ~ -350 
km s _1 ). The red wing (v > km s _1 ) of the feature extends 
up to v ~ 50 km s _1 , while the blue wing (v < km s -1 ) is 

much stronger and expands as far as v 410 to -730 km s . 

This makes the [Oi] 6300A line in these sources outliers in 
Figs.0and|S] Because of the obvious differences, the parame- 
ters describing the [O i] profiles of these three sources are not 
included in the computation of the average [Oi] parameters. 
The asymmetry parameter on the other hand is close to 1 for all 
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Fig. 10. The FWHM of the detected [O i] emission lines versus 
the asymmetry parameter of the feature. The plotting symbols 
are as in Fig. [7] The dashed horizontal line represents the mean 
value (1.04) fartt, the vertical line indicates the mean FWHM 
(47 km s _1 ). The typical error bar is represented in the upper 
right corner. 



three sources. Indeed, the features are quite symmetric around 
the (strongly blueshifted) centroid position. This is due to the 
fact that the blueshifted emission peak in these features is of 
comparable strength to the emission peak near km s" 1 . 

For 30 objects, the 6363A line region is covered by our 
spectra. In 15 sample sources, where the [Oi] 6300A line was 
detected, the [Oi] 6363A line was detected as well. The profile 
of the [Oi] line in these objects is pure-emission. Six sources 
for which the [Oi] 6300A line was detected, did not display 
the 6363 A line. In the remaining 9 objects, no [Oi] lines were 
detected. This includes all 5 sources with an absorption-line- 
rich 6300A wavelength region, for which we have data around 
6363A. The mean ratio of the EWs of the 6300A and 6363A 
lines EW(6300)/EW(6363) is 3.1+0.9, which is within the er- 
ror equal to the ratio of the Einstein transition rates of the lines 
(^6300/^6363 = 3.0). The lower limits derived for the sources 
with an undetected 6363A line are also consistent with this 
value. The mean difference between the FWHMs of both lines 
is close to km s _1 , FWHM(6300A)-FWHM(6363A) = 3+20 
km s _1 . The [Oi] 6363 A emission line has the same origin as 
the [Oi] 6300A line, as one would expect for two lines that 
originate from the same upper energy level ('Z^)- The spec- 
tral information in one line is exactly the same as in the other. 
The only parameter that can influence the EW ratio of both 
lines is differe ntial extinction. U sing a standard interstellar ex- 
tinction law dFluks et alJll994h and the error bar on the EW 
ratio (~30%), one can derive that the visual extinction towards 
the emission region must be less than 30. Also for Z CMa — 
in which the profile is very broad, and we have a sufficient 
S/N in both the 6300 and 6363A line to allow an accurate 
comparison — the two [O i] lines have comparable profiles. The 
[S n] emission line at 6731 A -which is present in the WHT 
spectrum of Z CMa- has a similar shape as the [O i] lines, be 
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Table 3. The different categories of Ha profiles versus the 
types of [Oi] profiles. The entries of the table are the num- 
ber of stars belonging to the corresponding category. U: unde- 
tected feature, E: pure-emission profile, A: 6300A region rich 
in photospheric absorption lines, S: single-peaked profile, D: 
double-peaked profile, P: P Cyg-like profile, inv. P: inversed 
P Cyg-like profile. 



[Oi] 

Ho- 


U 


E 


A 


Tot. 


U 





1 





1 


s 





5 


5 


10 


D 


6 


15 


4 


25 


P 


4 


7 


1 


12 


inv. P 





1 





1 


Tot. 


10 


29 


10 


49 



it with a much weaker blue wing. This forbidden line is unde- 
tected (< 0.01A) in the other FEROS and WHT spectra. 

The wavelength region around the forbidden oxygen line at 
5577 A is covered by the FEROS and WHT spectra. No features 
where detected, with upper limits of 0.01 A. Assuming that the 
emission of both the 5577A and 6300A line are thermal, one 
can estimate the theoretical Saha-Boltzmann 5577/6300 inten- 
sity ratio. In our sample we measure upper limits for this ra- 
tio of the order of ~l/20. The observations hence exclude that 
the oxygen lines are produced by thermal emission of oxygen 
atoms at temperatures above ~3000K. This is a strong indica- 
tion that the source of the [O i] emission at 6300A in HAEBEs 
cannot be found in thermal excitation of oxygen in a "super- 
heated" surface layer, as is the case for T Tauri stars (see later). 

The Ha features in our sample were classified based on the 
observed profile: single-peaked, double-peaked, P Cyg-like or 
inversed P Cyg-like. Despite the object's name, in our spectrum 
of LkHa 225 no Ha feature was detected. Since this spectral 
feature has previously been o bserved in emission in both com - 
ponents of this binary system (Maaakvan & Movsesvan 1997), 
the absence of the line in our spectrum suggests Ha variabil- 
ity. 10 of the 49 sample sources (20%) display single-peaked, 
25 (5 1 %) double-peaked and 12 (24%) P Cyg-like Ha emission 
profiles. In S V Cep, an inversed P Cyg-profile is observed. This 
is in good agreement with the H a distrib ution in the s ample s 
of iFinkenzeller & Mundll dl984 and iBohm & Catalal dl994l) : 
50% double-peaked, 25% single-peaked and 20% P Cyg-like. 
In Table [5] the different categories of the Ha profile are listed 
versus the types of the [Oi] profile. When plotting the [Oi] 
luminosities of the sample stars versus the Ha luminosities, 
regardless of the type of profile, a clear correlation is noted 
(Fig. [TTJ- Z CMa, PV Cep and V645 Cyg seem to have sig- 
nificantly more [Oi]-to-Ha luminosity than other sources. No 
significant differences between the other groups are noted. The 
average luminosity ratio L([Oi])/L(Ha) is 1.4 x 10~ 2 . 

3.2. [Oi] versus the SED parameters 

In this section we compare the parameters describing the [O i] 
6300A emission line to the SED parameters. In Fig.ll2lthe [O i] 
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Fig. 11. The [Oi]-versus-Ha luminosity. The sources with a 
negative Ha EW have been included. The typical error bar is 
presented in the upper left corner. Plotting symbols are as in 
Fig. The dotted line represents the median luminosity ratio 
L([Oi])/L(Ha) = 6.3x 10~ 3 . 
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Fig. 12. The [O i] 6300A luminosity versus the stellar UV lu- 
minosity. The plotting symbols are as in Fig.0 The typical er- 
ror bar is indicated in the upper left corner. The dotted line rep- 
resents the median luminosity ratio L([Oi])/L uv = 1.5 X 10~ 5 . 

luminosity is plotted versus the UV luminosity of the central 
star. The dashed line represents the median [O i]-over-UV lumi- 
nosity ratio (1.5 x 10~ 5 ) of the detected lines, while the average 
value is ~3.9 xl0~ 5 . In general, the [Oi] luminosity increases 
with increasing L uv . The different plotting symbols represent 
the three groups of HAEBE stars in our sample. Note that there 
is a considerable number of (almost all group II) sources with 
an undetected [O i] 6300A feature. The upper limits clearly de- 
viate from the average. 

The group III sources display the strongest [O i] emission. 
All 5 sources in this group have high L([Oi]) in the range be- 
tween 2.7x 10~ 3 and 3.8L . The median [O i] 6300A luminosity 
is 6.5 x 10 _1 L Q . Of the 12 group I objects in the plot, only one 
source does not show a detected [O i] 6300A emission line. The 
median [O i] luminosity of the detected features in this group is 
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5.9 x 1O~ 4 L . About 43% of the group II sources in this sam- 
ple (excluding those with absorption-polluted spectra) do not 
have a detected [O i] emission line. Moreover, even the median 
luminosity of the detected features in group II (3.2 x 1O _4 L ) 
is smaller than in the other two groups. Fig.ll3lshows the his- 
togram of the EWs of the group I and II sources. The group I 
sources have slightly stronger [Oi] intensities than group II 
sources. Also the typical width of the features differs between 
group I and II. The latter is illustrated in Fig.^] The group II 
profile is ~20 km s^broader than the group I profile. We con- 
clude that there is a clear correlation between the IR-SED clas- 
sification of the sources and the strength and shape of the [Oi] 
emission line. 




1 



1 



0.0 0.5 1.0 1.5 2.0 

EW [A] 

Fig. 13. Histogram of the EWs containing the 13 group I and 
24 group II members for which we have a measurement of, or 
an upper limit to the EW of the [O i] 6300A emission line. The 
filled part of the bars represents group I sources, the open part 
group II members. The bin size is 0.05A. The group I sources 
have slightly higher EW values than their group II counterparts. 



In Table@]we have summarized the parameters describing 
the [Oi] emission lines. The sample sources are listed accord- 
ing to their classification. 

3.3. [Oi] versus the infrared solid-state bands 

In this section we look for correlations between the PAH pa- 
rameters and the [O i] 6300A parameters for the 40 sample ob- 
jects which were included in the analysis of AV04. 

Fig. shows the estimated PAH luminosity versus the 
[Oi] 6300A luminosity. The 31 objects that are included in this 
plot have a pure-emission [O i] profile. A clear positive corre- 
lation is seen for the 14 stars that have detected PAH as well 
as [Oi] 6300A emission. The PAH luminosity is on average 
350 times larger than the [O i] luminosity. All group III sources 
in this sample have strong [O i] luminosity, but weak or unde- 
tected PAH emission. The PAH and [Oi] luminosities of the 
group II sources are weak or undetected in most cases, while a 
typical group I source displays strong PAH and [O i] emission. 
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Fig. 14. The average [O i] 6300A profiles for group I and 
group II. The horizontal dotted lines are plotted at and 0.5, 
while the vertical dotted line represents the centroid position 
of the feature. The 12 group I and 12 group II sources which 
display [O i] emission are included. The profile of each source 
was shifted to its centroid position, continuum-subtracted and 
normalized by dividing by its peak-flux. The group II profile 
(dashed line) is ~20 km s 'broader than the group I profile (full 
line). 
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Fig. 15. The PAH luminosity versus the [O i] 6300A luminos- 
ity. The plot includes the 12 group 1, 14 group II and 5 group III 
sources for which (an upper limit of) both measurements was 
available. The dotted line represents the median PAH-over-[0 1] 
luminosity ratio L(PAH)/L([Oi]) = 130. Plotting symbols are 
as in Fig.0 The typical error bar is indicated in the upper left 
corner. 



When comparing the luminosity of the amorphous 10 mi- 
cron silicate feature to the [Oi] 6300A luminosity (Fig. I16> . 
no clear correlation is seen. The arrows which indicate upper 
limits support this non-correlation. 
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Table 4. The [Oi] parameters for each sample source. The objects are listed according to their IR-SED classification. The 
distance to the source, the effective temperature, UV luminosity and bolometric (= integrated observed) luminosity are giv en. 
References for the firs t two values can be found in Table 6 of AV0 4, unless otherw ise indicated: A de Zeeuwetal. (199j; B 
lFinkenzellerl(ll985h: c lGrav & Corballvl(ll993l) : D lGrav & Corballvl(ll998l) : £ lHoukl(ll978l) : F lMora et all feOOll) : G IVieira et alJ 
(2003); H van den Ancker et al. ( 1998). The type of [Oi] 6300A emission profile observed in each object has been indicated; 
E: pure-emission profile, A: 6300A region rich in photospheric absorption lines, U: undetected feature (upper limit), n.s.: no 
spectrum available. The typical errors on the FWHM and centroid position are ±8 and ±15 km s -1 respectively. The thirth column 
contains the stellar radial velocity em pl oyed in this study. This value is es tim ated from the available spectra, unl ess otherwise 
indicate d: a lArce & Goodmanl J2002I): b lArellano Ferro & Giridharl J2003I): c barbier-Brossat & Figonl d2000l) : d bunkin et all 
Jl997al) : e IPietu et al.1 d2003l) : / lReipurth et alJ dl996l) : g lYonekura et alJ Jl997l based on CO lines). f Based on Ca n fines. O 
Average of radial velocities of stars in the same field (radius~2')- The last two columns of the table give the type and EW of the 
observed Ho- profile; S: single-peaked; D: doubl e-peaked; P: P Cygni-lik e; inv. P: inversed P Cyg-li ke. When no spectrum w as 
available, literature valu es were used: /] Arellano Ferro & Giridharl ll2003l): ' Ide Winter & Thel Jl990l) : ^ bunkin et all dl997ah : k 
|PunkinetalJlll997bll : ' iFinkenzeller & Mundll i ll 9841) : m lMerin et alJll2004 : " IPerez & Gradvl ill 9971) . 
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4. Interpretation 



A few sources like the group III outliers Z CMa, PV Cep and 

xts' A e 1 i. t t irv ^nnnr _i tt 



14 



Acke, van den Ancker & Dullemond: [O i] emission in HAEBE disks. 



10' 




1 1 


1 


1 


10 1 

>^ 




♦ 




T 

•is 




| 10° 

Z3 
0) 

B 10_1 




o * o^- 
f o *♦ 


♦ 


I 

I 


o 

CD 










£ io~ 2 




' I 






O 

^~ 10~ 3 




,i ,4- i 


1 





10 5 10 4 10 3 10 2 10 1 10° 



[01] luminosity [L sun ] 

Fig. 16. The luminosity of the amorphous 10 micron silicate 
feature versus the [O i] 6300A luminosity. The plot includes the 
12 group I and 15 group II sources for which (an upper limit 
of) both measurements was available. The group III sources are 
left out since they display silicate absorption at 10/um. Plotting 
symbols are as in Fig. The typical error bar is indicated in 
the lower right corner. 
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Fig. 17. The [O i] 6300A emission line in the spectra of HD 
100546. The double-peaked profile is clearly seen in the ESO 
3.6m and CAT spectrum. Due to the lower spectral resolution, 
the profile remains unresolved in the FEROS spectrum. The 
horizontal dotted lines represent the continuum levels of the 
three spectra. The vertical line is the mean centroid position of 
the feature at 3 km s _1 . 



source HD 141569 have broad features (FWHM >100 km s _1 ) 
with pronounced blue wings. A few more sources have tow- 
velocity blueshifted centroids. Nevertheless, the majority of 
observed [O i] 6300A emission lines have narrow (FWHM ~ 50 
km s _1 ) symmetric profiles, centered around the stellar radial 
velocity. Some of the high-resolution spectra display double- 
peaked profiles with a peak-to-peak distance of ~10 km s . 
The low velocities, symmetry of the feature and the peak-to- 
peak separation correspond to what one would expect from an 
emission-line region at the circumstellar-disk surface. We in- 
terpret the observed [O i] lines in the majority of our sample as 
being circumstellar-disk emission features. The forbidden-line 
emission region is located in the disk's atmosphere; a warm 
layer, directly irradiated by the central star and corotating with 
the disk. The blue wing which is present in a minority of the 
cases is most probably emanating from an outflow, with the red 
wing eclipsed by the circumstellar disk. In this section we will 
discuss the arguments for this interpretation. 

When investigating the narrow profiles of the non- 
blueshifted features in the high-resolution ESO 3.6m spectra, 
a double-peaked profile is seen in a some cases. In Fig. El 
this specific shape in the spectrum of HD 100546 (group I) 
is shown. The two peaks are located at equal distance (~ 
+6 km s _1 ) from the centroid position. Simila r profiles have 
been observed in CO emission lines (e.g. Cha ndler et alll993t 
iThi et all200l and have been attributed to the Keplerian rota- 
tion of the circumstellar gas disk. Also in group I sources HD 
97048, HD 135344 and T CrA, for which we have ESO spec- 
tra, this double-peaked shape is clearly seen. The other spectra 
have a lower spectral resolution, and thus the double peaks are 
less clear. Considering the latter, more sources with symmet- 
rical low-velocity profiles, like AB Aur, might have a double- 



peaked line shape. In section [5] we model the spectral profiles 
of some of these features. 

A second indication that the forbidden-line emission region 
is connected to the circumstellar disk lies in the correlation 
between the strength of the [Oi] line and the SED classifica- 
tion of the sources. As explained before, group I sources are 
suggested to have a flared-disk geometry, while it is believed 
that group II members have self-shadowed disks. The disk ge- 
ometry can play an important role in the amount of neutral 
oxygen in the upper level of the 6300A line. Apart from ther- 
mal excitation, photodissociation of OH and H2O molecules by 
UV photons can significantl y increase the num ber of excited 
oxygen atoms Jvan Dishoeck & Dalgarnd| 19841) . When these 
molecules are exposed to high UV fluxes, e.g. at the surface of 
a flared disk, the non-thermal population of the upper level of 
the 6300A line occurs very efficiently. The outer parts of self- 
shadowed disks lie completely in the shade of the puffed-up 
inner rim. Consequently, the UV radiation of the central star 
cannot reach the disk surface and no significant photodissocia- 
tion of OH and H2O molecules will take place. This would re- 
sult in weaker [0 1] emission for group II sources. The absence 
of the feature in 40% of the group II sources and Fig.ll3lindeed 
show that the correlation between [0 1] 6300A luminosity and 
M01 classification exists, as described in SI3.2I 

Furthermore, the [0 1] luminosity is also correlated to the 
PAH luminosity of the source (see Fig. I15> . These molecules 
need stellar UV photons to get excited. In flared disks, signif- 
icantly more PAH emission is observed than in self-shadowed 
disks (AV04). The observed correlation between the [Oi] and 
PAH luminosity suggests that both the oxygen atoms and the 
PAH molecules radiate from the same location: the disk's atmo- 
sphere. This low-density, UV-photon-immersed region at the 
disk's surface offers the ideal locus for the photodissociation of 
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OH and H2O and the consequent 6300A forbidden-line emis- 
sion. 

The non-correlation of the [Oi] 6300A luminosity with 
the luminosity of the amorphous 10 micron band is also ex- 
pected under the present hypothesis. The small silicate grains 
that cause this emission feature are thermally excited and do 
not need direct stellar flux. The observed lack of correlation is 
similar to the non-correlation between the PAH luminosity and 
the luminosity of the 10 micron feature as described by AV04. 




5. Modeling the [Oi] emission region 

In this section we model the [0 1] emission in the atmosphere 
of a flared disk. We have co mputed the structure of a flared 
disk JChiane & Goldreichl 19971) and determined the layer from 
which the [Oi] emission emanates. In this layer, the density 
is low and the UV flux abundant. We calculate the intensity 
of the [Oi] emission for different stellar and disk parameters. 
In order to model the line profiles observed in the spectra, the 
theoretical emission-line profile was determined by convolving 
the computed intensities with a Keplerian-rotation profile. 



5.1. The flared-disk model 

We have implemented the flared-disk model of 
I — I J — ' — 1L 

Chiang & Goldreich ( 1997) with some improvements de- 
scribed by Dullemond et alJ 11200 ll their section 2.1.1 and 
2.1.2). The input values of the model consist of the stellar 
parameters (M+, and L*), dust opacities k v and the surface 
density S. The latter quantity is assumed to be a power-law 
function of the radius R to the star: Z = E CKfAU])"^. The 
surface density at 1 AU (Sq) and the power (J3) can be chosen 
freely. The structure of the disk's interior is then calculated 
iteratively by demanding vertical hydrostatic equilibrium at 
each radius. The outp ut quantities (in the nomenclature of 
iDullemond et alJ EoOl ) include the pressure scale height h cg , 
the disk surface height H cg , the midplane temperature 7 1 , and 
the surface temperature T s . These are all a function of the 
radius R. The density p is a Gaussian function in the vertical 
(i.e. z) direction, centered around the midplane (z = 0). Its 
half-width depends on the pressure scale height h cg , which is a 
function of the midplane temperature T t . 

The temperature distributio n in the disk is cal culated 
more precise than described in IChiang & Goldreichl ill 9971) 
and IDullemond et alJ J200ll) . These authors assume a two- 
temperature model which is determined by the midplane 
temperature T, and the temperature at the optically thin 
surface layer T s . Inspired by t he full-fletc h ed, co mpu- 
tationally demanding model s of IDullem ond (2002) and 
IDullemond & Dominikl J2004I) . we allowed for a temperature 
gradient in vertical direction. The stellar flux penetrates the 
disk, directly heating the circumstellar matter. At a certain ra- 
dius R from the star and height z above the midplane, the tem- 
perature T(R, z) can be determined using 



R [ALU 



Fig. 18. The temperature structure of a flared-disk model. On 
the x-axis the radial distance R to the star is plotted, on the y- 
axis the ratio of the height above the midplane z over R. The 
color scale is linear and ranges from to 3500K. The dashed 
line represents the disk surface height H cg . 

r^ r°° 

J fi v (r*)exp(-T rad OR,z))K v £/v (2) 

in which and are the stellar radius and temperature. The 
radial optical depth r ra d can be approximated by the ratio of the 
vertical optical depth r vert and the flaring angle of the disk a: 



TrzdiR, Z) = 



a(R) 



(3) 



Jo 



B v (T(R,z))Kydv = 



This approach is similar to the one used to compute the sur- 
face scale height of the disk in the original lchiang & Goldreichl 
model. At the surface, where T ra( j <k 1, the temperature T(R,z) 
is equal to the optically thin T S (R). Deeper in the disk's in- 
terior, the T(R,z) computed from equation (0 drops until it 
reaches the midplane temperature Tj(R). At these locations, the 
dominant heating source is not the direct stellar flux, but the 
flux radiated downward from the disk surface. The latter in- 
direct flux maintains the midplane temperature Tj and is only 
important in regions were no direct stellar flux can penetrate. 
Therefore, we have set all T(R,z) equal to T,(R) at these loca- 
tions. In Fig.[D5]the temperature structure of a typical flared- 
disk model is shown. 

We note that the temperature stratification in the disk and 
the density function are not computed self-consistently. The 
temperature gradient is calculated a posteriori; the effect of 
the higher temperatures in the upper layers of the disk is not 
included in the determination of p. The deviation from the self- 
consistent calculation is however very small. 

5.2. The [Oi] emission 

To determine the [Oi] 6300A emission emanating from a 
flared circumstellar di s k, we have implemented the model of 
LStorzer & Hollenbachl (2000) for [Oi] emission. They have 
modelled the optical forbidden-line emission from a plane- 
parallel semi-infinite photodissociation region (PDR). Both 
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thermal and non-thermal emission are included. Slight changes 
were applied to the method to adapt it to our specific case. 

5.2.1. Thermal [Oi] emission 

Following IStorzer & Hollenbachl fcOOOl) . a five-level oxygen 
atom is assumed. The atomic data and references can be found 
in their paper (Table 1 of the Appendix). Since we are only in- 
terested in modeling the 6300A line, the relevant transitions 



for the present paper are 
3 P 2 (44.2 //m), >D 2 



J Pj — -> 5 P 2 (63.2 jum), J P — > 
3 P 2 (6300.3A) and l S l D 2 
(5577. 4A). Collisions with free electron s and atomic hydrogen 
are co nsidered. From PDR models (e.g. Hollenbach & Tielens 
1999, and references therein), one derives that the transition 
region from the atomic-H dominated upper layers to the disk 
interior where H 2 is abundant occurs at about Ay - 1. Based 
on the latter, we assume that only at Ay < 1 free electrons and 
H atoms are important collisional partners for the O 1 atom. In 
optically thick regions, the collisional rate is set to zero. 

The proton number density n p depends on the radial and 
vertical position (R, z) in the disk and the input surface den- 
sity S. It is computed from the flared-disk model. The oxygen 
number density is the product of the proton density and the 
fractional oxygen abundance 



n(O) = n p 



(4) 



A typ ical interstellar value for [O/H] is 5 x 1(T 5 (e.g.. lAlleiJ 
|l976J). In the region where oxygen is not ionized, the atomic 
hydrogen is neutral as well, because the ionization energy of 
the two species is comparable. In this region the free electrons 
are almost all due to the ionization of neutral carbon atoms. 
Assuming all carbon atoms in the considered region are singly 
ionized, the free-electron number density is 



(5) 



A typical interstellar fractional abundance of carbon is 3x 1CT 4 . 

To determine the thermal [0 1] 6300A emission, the popu- 
lation of the upper level of this line ( l D 2 ) needs to be known. 
Using the formula of thermal equilibrium (see Appendix), one 
can determine the relative population of two levels a and b for a 
temperature T. From the latter ratios, the relative fraction F of 
O 1 atoms in the l D 2 state can be calculated. The number den- 
sity of the thermally populated upper level of the 6300A line is 
n u ,ih = Fn{0). 

Since the typical temperature range (10-1500K) in our disk 
models is rather low to thermally excite the oxygen atoms, the 
emanating [0 1] 6300A emission is very weak. The intensity of 
the thermal emission is more than a million times smaller than 
the non-thermal [0 1] emission, which is discussed in the next 
section. 

5.2.2. Non-thermal [Oi] emission 

Following IStorzer & Hollenbachl J2000h . the dominant non- 
thermal excitation mechanism for neutral oxygen atoms is the 



photodissociation of OH molecules. This results in a hydro- 
gen atom and an excited oxygen atom. A fraction of the lat- 
ter (~55%. lvan Dishoeck & Dalgarno 1984) find themselves in 
the upper state ( l D 2 ) of the 6300 A line. This mechanism hence 
produces strong non-thermal [Oi] emission in regions where 
the photodissociating UV flux is abundant and the densities 
high enough to have a sufficient amount of emitting oxygen 
atoms. 

In our simple model for non-thermal emission, we assume 
that the fractional OH abundance e{OH) is constant through- 
out the disk. In reality, this may not be the case, as the opti- 
cally thick disk interior will have a much higher e(OH) than 
the photon-immersed surface layers. Nevertheless this assump- 
tion needs not to be valid for the entire disk, but only for the 
[Oi] emission region. The latter will be located close to the 
t = 1 surface and is geometrically quite thin, because the 
proton density drops off rapidly when moving away from the 
disk midplane while the high optical depth in the disk inte- 
rior prohibits photodissociation of the OH molecules. From 
the OH number density n(OH) = e{OH) x n p , the den- 
sity of non-thermally-excited oxygen atoms can be determined 
dStorzer & Hollenbachl 19981 adapted to the UMIST rate coef- 
ficients): 



n, 



u,non—th 



/Go 



3.50 x 1(T 



n(OH). 



(6) 



7 In this formula, / (« 55%) is the fraction of oxygen atoms 
released by the photodissociation of OH in the upper level of 
the 6300A line. The sum (A) of the Einstein transition rates 
from the 1 Do-leve l downwards is approximately 8.4 X 10 3 s _1 
dOsterbrockll989l) . 

The local [Oi] 6300A emissivity j(R,z) (in erg s _1 cm -3 
sr" 1 ) is given by 



j(R,z) = ^-hvA b300 n u (R,z) 
An 



(7) 



The E instein coefficient A^oo is 6.3 x 10~ 3 s _1 ( Oster brockl 
1989). The value hv is the energy of a 6300A photon. The den- 
sity n u is the sum of the densities of the thermally and non- 
thermally excited oxygen atoms. In Fig.[H|] j(R,z) is plotted 
as a function of the vertical and radial location in the disk. The 
dashed line represents the disk surface height H cg of the disk. 
The emission is located around this region. 

Integrating the emissivity over the vertical direction, con- 
sidering the optical depth T6300 at 6300A, one extracts the sur- 
face intensity I(R) at each radius R: 



I(R) 



-I 



j(R, z) e~ T6m dz 



(8) 



The unit of I(R) is erg s" 1 crrT 2 sr -1 . Note that, due to the shal- 
low angle under which the stellar light impinges on the disk, 
one looks much deeper into the disk in the vertical than in the 
radial direction. When looking at the disk face-on, the observer 
practically sees the entire [0 1] emitting region (on one side of 
the disk). 
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Fig. 19. The [O i] 6300A emissivity in the flared-disk model. 
The x-axis gives the radial distance R to the central star in AU, 
the y-axis plots z/R, the ratio of the vertical height and the ra- 
dial distance. The colors represent the local [Oi] 6300A emis- 
sivity on a log-scale which covers 10 orders of magnitude. On 
the color bar, the emissivity is given in [erg s _1 cirT 3 sr~']. The 
disk region sampled by the model ranges from roughly 0.6 to 
1.5 times the disk surface height H cg . The latter is indicated by 
the dashed line. The geometrical full width at half maximum of 
the [Oi] emission region in the vertical direction is 5-10% of 
H cg . 

5.3. The line profiles 

To convert the intensity-versus-radius function I(R) to an 
intensity-versus-velocity profile, we assume that the rotation 
of the disk is Keplerian 3 . We anticipate this to be a fairly ac- 
curate assumption, since the disk mass is ex pected to be much 
smaller than the mass of the central star (e.g., Ac ke et al .2004). 
Additionally we adopt that the Keplerian rotation velocity of 
matter in the upper layers of a flared disk is the same as the 
rotational velocity in the midplane. Fig. [^illustrates how the 
profile is computed. The light-grey band in the pictogram rep- 
resents a circular orbit in the inclined -geometrically flat- disk 
where the [Oi] intensity is constant. The disk's matter in this 
band rotates around the central star with a Keplerian velocity 
Vkep = VGM*//? or bit in which R OI bk is the radius of the orbit 
and M* is the mass of the central star. The narrow darker- 
grey band represents the parts of the disk that move towards 
(or away from) the observer with the same projected velocity 
Vprqj • This area is the region where 

GM+ , , 

R = sin 2 i sin 2 9 (9) 

Vproj 

with 9 the angle in the disk plane (9 - 0° towards the observer). 
The inclination i is defined to be 0° for pole-on disks. The black 

3 The IDL code keprot.pw, which converts an intensity-versus- 
radius profile into an intensity-versus-velocity profile, can be 
downloaded at 

http://www.ster.kuleuven.ac.be/~bram/OI/keprot.pro 
http://www.ster.kuleuven.ac.be/~bram/OI/keprot.README 




projected velocity 



Fig. 20. Pictogram of a flat rotating circumstellar disk, seen 
under an inclination of i — 45°. The location of the central star 
is indicated by the asterisk. The light-grey band represents a 
circular orbit. In this orbit, the rotational velocity as well as the 
[Oi] intensity (I ) are constant due to the assumption of ax- 
isymmetry. The narrow dark-grey band represents the region 
where the projected rotational velocity is equal to a certain 
value v p . At the cross-section of both surfaces the projected 
velocity is equal to v p and the intensity is equal to I e . The in- 
sertion shows the intensity-versus-velocity profile of the light- 
grey band. The dark-grey part represents the contribution of the 
cross section of the dark-grey band and the orbit in the disk. 

surface is the cross-section of both bands. By adding up these 
regions and multiplying with the intensity, the intensity-versus- 
velocity profile is built up for each projected (=observed) ve- 
locity. 

The computed theoretical profile is then convolved with a 
Gaussian function 4 . This function mimics the effect of instru- 
mental broadening. The width of the -unresolved- telluric lines 
in each type of spectrum were used as the width of the instru- 
mental Gaussian profile. 

5.4. Discussion of the model parameters 

The input parameters for our model that calculates the 
non-thermal emission consist of some stellar parameters 
(L*, L uv , M^, T*) which are derived from the photometric data. 
Furthermore, a dust opacity tabl e and dust-to-gas mass ratio 
are needed. Following the bullemond etail £2001) model for 
AB Aur, we assume that the dust consists of olivines and rep- 

4 The IDL code convolve. pro, which convolves an intensity-versus- 
velocity profile with a Gaussian function, can be downloaded at 
http://www.ster.kuleuven.ac.be/~bram/OI/convolve.pro 
http://www.ster.kuleuven.ac.be/~bram/OI/convolve.README 
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Fig. 21. The intensity-versus-radius distribution I(R) for four 
template models, with different values for the surface-density 
parameters. From top to bottom, /3 is 1.0, 2.0, 2.5 and 3.0 re- 
spectively. The value of So is chosen such that the disk masses 
of the four models are the same. The intensity function is mul- 
tiplied by the total surface S (R) of the ring at radius R. In this 
way the fluxes from the inner and outer parts of the disk can 
be compared directly. The vertical dashed lines mark the radii 
where the Keplerian velocity is 5, 10, 50 and 100 km s . The 
shoulder of the intensity distribution shifts closer to the star 
with increasing /3. Furthermore, the total emission is smaller 
when the power-law describing S has a steeper slope. 



resents 1% of the total disk mass. The stellar parameters are 
Li, = 47 L , L uv = 36 L e , M* = 2.4 M G and = 9520K. 
The free parameters in our [Oi] emission model are the sur- 
face density S = So (R[AU])~P, the fractional OH abundance 
e(OH) and the disk's inclination i. In the template model, we 
take S = 10 4 g cnr 2 , f3 = 2, e(OH) = 1 x 10~ 6 and i = 45°. 
The inner and outer radius of the disk, /?,„ and R out , are 0.1 
AU and 100 AU respectively. We discuss the influence of the 
four free input parameters on the shape of the emission profile, 
starting with the surface density. 

The two parameters So and /3 describing the power-law sur- 
face density are coupled when the disk mass is known: 



M, 



disk 



= 2n f " 



So R l ^dR 



(10) 



The disk mass of our template model is M c u s k = O.O49M . 
In Fig. the intensity-versus-radius distribution of the [Oi] 
6300A emission is plotted for four models. The input param- 
eters of the latter are the template values, except for S. The 
power /3 ranges from 1.0 to 3.0, while So is appropriately 
adapted to ensure that the disk mass remains unaltered. The 
total emission intensity decreases with increasing /3. Fig. [22] 
shows the line profiles, corresponding to the intensity distribu- 
tions in Fig. 1211 

The fractional OH abundance e(OH) scales linearly with 
the total intensity of the profile. It does not alter the shape, as 
we assume it is constant throughout the [O i] emission region. 

Observations of Doppler-broadened spectral lines do not 
observe the real, but the projected velocities in the system. In 
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Fig. 22. The line profiles /(v) corresponding to the intensity- 
versus-radius distribution in Fig. |2^ From top to bottom, /3 is 
1.0, 2.0, 2.5 and 3.0 respectively. Since the intensity-versus- 
radius distribution falls off rapidly in the models with large [} 
values, the corresponding profiles have peaks at higher veloci- 
ties than the profiles belonging to low-/? models. 
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Fig. 23. The line profiles I(v) of the template model seen un- 
der different inclinations. From top to bottom, the inclination 
is 18°, 30°, 45° and 70° respectively. Note the decrease in inte- 
grated intensity and the shift of the peaks with increasing i. 



the present flared-disk model, the inclination i is a free param- 
eter which alters the shape as well as the intensity of the com- 
puted profile. The total integrated intensity J I(v)dv is propor- 
tional to cos i, while the position of the two typical peaks of 
the line changes with sin/. Fig. 1231 shows the line profiles of 
the template model as seen under different inclinations. Note 
that inclinations larger than about 70° are not relevant in flared 
disks, since at such high i, the outer parts of the disk would 
occult most of the disk surface. 

6. Comparison with the observations 

In this section we compare our model results to the observa- 
tions. We focus on the sample stars which display a narrow, 
symmetric and centered profile. Since the thermal emission is 
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Fig. 24. Similar plot as Fig.^] The filled diamonds are the ob- 
served [O i] intensities for the group I sources. The open circles 
represent model results. The typical error bar for the observa- 
tions is given in the upper left corner. The models differ from 
each other in stellar parameters; from the lower left to the up- 
per right, the model [Oi]-versus-UV luminosity for a typical 
FOV, A3V, AOV, B8V, B5V and BOV star are plotted respec- 
tively. The error bars on the model determinations represent 
the spread in results when simultaneously varying the surface 
density £ and inclination i in the ranges described in Figs. |22] 
and[23] The dotted lines represent the model results when alter- 
ing the fractional OH abundance e(OH) and using the average 
[O i]-versus-UV luminosity ratio of the models. The abundance 
is indicated next to each line, with a{-b) representing a x 10~ & . 
The fractional OH abundance of the models is 1 x 10~ 6 . 



weak compared to the non-thermal emission for the relatively 
cool disks studied in this paper, we will only consider the non- 
thermal emission component in the following discussion. 

We compare the observed [O i] 6300A emission luminos- 
ity and its theoretical counterpart to the UV luminosity in 
Fig. |24] Our model applies for fiared-disk geometries, there- 
fore only the group I sources have been plotted in the figure. 
The model results and the observations agree nicely for an ac- 
ceptable range of the free parameters X, i and e(OH). Note that 
all models in the figure have the same disk mass as the template 
model. Varying the disk mass affects the [Oi] 6300A emission 
in a comparable manner as when the surface density is altered. 
Fig-EJcompares the [Oi] luminosity to the effective tempera- 
ture of the group I objects and the models. Again, reasonable 
values for the free parameters lead to theoretical [O i] luminosi- 
ties close to the observed values. 

Not only the strength of the [O i] 6300A emission can be re- 
produced with our model. After scaling the total theoretical to 
the observed intensity, the typical FWHM and double-peaked 
shape of the narrow centered features can be explained by as- 
suming Keplerian rotation of the flared disk. In Fig.[26]we have 
overplotted observed profiles with theoretical ones. For each 
star, we have used the stellar parameters (L*, T*, M*) as input 
parameters. The remaining parameters were set to the template 
values, unless otherwise indicated in the plots. Note that the 
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Fig. 25. The [O i] 6300A luminosity versus the effective tem- 
perature of the group I sources (filled diamonds). The open cir- 
cles and dotted lines plotted in this figure represent the same 
models as in Fig. 1241 

theoretical profiles are not fitted to the data, but are just over- 
plotted model profiles, rebinned to the spectral resolution of 
the observed spectrum. The similarity between theory and ob- 
servation is nevertheless striking. For AB Aur and HD 100546, 
the in clinatio n of the disk h ave been determine d to be i - 76° 
(Mannings & Sarg ent! 19971) and i = 51° jAugereau et alJ200ll) 
respectively. The observed profiles of the [Oi] 6300A line can 
be approximated by theoretical line profiles for these inclina- 
tions when using e(OH) « 2 x 10~ 6 and yS = 1.5 for AB Aur 
and e(OH) « 6 x 10~ 7 and p = 2.5 for HD 100546. 

The flaring of the disk is due to the dust opacities, while 
the gas opacities are negligible. Since dust evaporates above 
the dust sublimation temperature T s m ~ 1500K, we have 
shifted the inner edge of the flared-disk from R m =0.1 AU 
to the radius where the surface temperature T s equals T m bi. 
For a typical FOV, AOV and BOV star, the inner radius hence 
becomes 0.2, 0.8 and 30 AU respectively. As a consequence, 
the high-velocity (|v| > 50 km s _1 ) wings of the theoretical 
[O i] emission line profile disappear, in agreement with the ob- 
served spectra of narrow 6300A features which do not display 
extended wings (e.g., AB Aur, HD 97048, HD 100546, HD 
135344, HD 169142, R Cr ATCrA, HP 1792 18, HD 190073). 
Furthermore, the models of lDullemondl J2002h suggest that the 
inner rim of the disk is puffed-up and casts its shadow of the 
first few AU in a flared geometry. This effect could be mim- 
icked in our flared-disk model by shifting the inner radius /?,„ 
even further out (up to ~5 AU). The effect of the latter on the 
theoretical line profile would be mostly a decrease of the wings, 
without reducing the total intensity much. 

The innermost parts of the circumstellar disk, where no 
dust can survive due to the high temperatures, are relatively 
small in spatial dimensions. Nevertheless, we cannot exclude 
that the contribution of these parts to the total [O i] emission 
strength is significant. Specifically, a few group II sources (e.g., 
VX Cas, V586 Ori, HD 95881, HD 98922, HD 101412, WW 
Vul, S V Cep) display broader, but relatively weaker [O i] pro- 
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Fig. 26. The observed spectra of group I sources AB Aur, HD 97048, HD 100546 and T CrA overplotted with the model results. 
The parameters of the model are adapted to the stellar parameters of each source combined with the template parameters, unless 
otherwise indicated in the upper right corner. We have shifted the theoretical profile to the centroid position of the observed 
feature. For HD 100546 two models with different parameter sets have been displayed to illustrate the difficulty to distinguish 
between models. The theoretical profiles extracted from our fiared-disk model fit the observed features strikingly well for a range 
of reasonable model parameters. a(-b) represents a x 10~ h . 



files. In our current interpretation and with no external bright 
UV source near, this [Oi] emission cannot emanate from the 
disk surface for the self-shadowing prohibits direct stellar flux 
to reach this area. In the latter cases, the emission might come 
from a rotating gaseous disk inside the dust-sublimation radius. 
The modeling of this inner gaseous region is however beyond 
the scope of the presen t paper. We refer to the recent paper of 
iMuzerolle et alJ l|2()04) for a more elaborate discussion on this 
subject. 

7. Conclusions and discussion 

The observed [O i] 6300A emission line in many HAEBE stars 
in our sample shows evidence for a rotating forbidden-line 
emission region. In this paper we suggest that the surface of the 
flared circumstellar disk around the group I objects is the per- 
fect location to harbor this emission region. The combination 
of the direct stellar UV flux and the relatively high densities 



in this region give rise to strong non-thermal [Oi] emission, 
which can explain the observed luminosities reasonably well. 
Furthermore, the shape of the spectrally resolved 6300A pro- 
file in the observations and the profile produced by our simple 
model indicate that Keplerian rotation indeed is the broadening 
mechanism for (at least the narrow component of) this line. The 
observed 6300A spectra of the group I members in our sample 
can be reproduced strikingly well. 

In the group of self-shadowed-disk sources, significantly 
more targets do not display the [Oi] 6300A emission line 
(43% versus 8% in group I). The line profiles of the detected 
[Oi] 6300A emission feature in group II are twice as weak 
as, and somewhat broader than the lines observed in group I 
spectra. In our current interpretation, the surface of a self- 
shadowed disk's outer parts is not directly irradiated by the cen- 
tral star. However, the non-thermal [O i] line formation mecha- 
nism — photodissociation of OH and H2O — may produce the 
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observed group II emission lines as well. Assuming that the 
[O i] emission emanates from the inner gaseous disk naturally 
explaines the weaker emission and the higher-velocity wings 
of the feature: the inner gaseous disk provides a smaller emis- 
sion volume and is located closer to the central star, where the 
Keplerian velocities are larger. 

In T Tauri stars, the less massive counterparts of HAEBE 
stars, the obser ved \Q il emission profile c an b e explained us- 
ing th e model of K wan & Tademarul (fr988) (e.g. lHartigan et all 
1 1995b . However, this model assumes that a "super-heated" disk 
atmosphere, fed by accretion, is present, in which the temper- 
atures are sufficiently high (~10 4 K) to produce thermal [Oi] 
emission. The absence of a strong UV excess and strong photo- 
spheric veiling, the relatively weak IR recombination lines and 
CO emission, and the presence of 10 micron silicate emission 
(AV04) in the majority of the group I and II sources in our sam- 
ple all indicate that these disks are passive. Furthermore, the 
observed upper limits for the 5577/6300 [O i] ratio indicate that 
the observed emission cannot be thermal. Each model which 
models the [Oi] emission based on thermal processes — in- 
cluding the disk wind model — will have difficulties reproduc- 
ing the strength of the [Oi] 6300A line in HAEBE stars. The 
non-detections of the [S n] 6731 A line in the spectra of group I 
and II sources in our sample confirm this picture. Moreover, 
amongst the late-type sample sources no strong [O i] emitters 
are present. In the T Tauri model the forbidden line emission 
is dependent on accretion rate and thus the model cannot ex- 
plain the absence of strong late-type [Oi] sources. The OH 
photodissociation model suggested in the present paper natu- 
rally predicts that passive disks around weak UV sources will 
not produce strong [O i] 6300A emission. 

The high-velocity blue wings in the [O i] 6300A line of a 
small minority of our sample stars cannot be accounted for 
by an emitting passive Keplerian disk. This emission feature 
is suggested to emanate from an outflow, of which the red- 
shifted part is occulted by the circumstellar disk. Note how- 
ever that this pronounced blue wing is often accompanied by a 
symmetric peak at low velocities (e.g., Z CMa, PV Cep, V645 
Cyg, HD 200775). The latter might again be formed in the 
surface layers of the rot ating disk. Alternatively, the model of 
iKwan & Tademarrl Jl98 8) for T Tauri stars may be valid for 
these objects: since the assumption of passivity is likely not to 
be valid for the disks of the group III members, these sources 
might resemble classical accreting T Tauri stars. As noted be- 
fore, accretion is needed to create the right settings for the T 
Tauri model to work. This idea is supported by the detection of 
the [Sn] 6731 A line in Z CMa, which is exclusive to this tar- 
get in the present sample. We note however that this behaviour 
seems to be rather exceptional for HAEBE stars. 

The values for the fractional OH abundance e(OH) needed 
to explain the observed [Oi] 6300 A luminosities are ~10~ 7 - 
10~ 6 . Observations of diffuse interstellar clouds show that 
relative OH abundances of this magnitu de (~1 0~ 7 ) occur 
in the interstellar medium (e.g., ICrutcherl 1979, and refer- 
ences therein). Nevertheless, these values are two orders larger 
than the abundances computed in recentjnodels inclu ding a 
full treatment of disk che mistry (e.g., iMarkwick et alj r2002: 
iKamp & Dullemondl2004l) . A possible reason for this discrep- 



ancy may be the exact location of the [Oi] emission region in 
our models: a shift to higher densities (i.e. to a lower vertical 
height z) would reduce the fractional OH abundance required to 
fit the observed [O i] intensity, since the latter is inversely pro- 
portional to the first. This effect can potentially be induced by 
the input dust opacities, which define the disk's flaring, but also 
the UMIST coefficients determine the exact location of the [O i] 
emitting region through formula l|6}. Alternatively, the chemi- 
cal models may be wrong due to uncertainties in the reaction 
rate coefficients or the incompleteness of the network. 

Boh m & Hirthl i 19971) have stated that t he forbidden-line 
emission region in Hillenbrand et group I sources 

cannot be located at the surface of the circumstellar disk, be- 
cause the suggested accretion disk would cover most of the 
forbidden-line emission region. The authors suggest that, even 
when the stellar light is able to reach the disk surface, a prob- 
lem remains: the upper layer would be geometrically very thin 
and the outer radius of the disk would need to be much larger 
than 100 AU in order to create enough emission volume to ex- 
plain the observed [Oi] intensities. Both problems are coun- 
tered when consid ering a flared disk with a puffed-up inner rim. 
Hillenbrand et al. have invoked a circumstellar accretion disk 
model to explain the observed near-IR (~2 pm) bump which is 
typical for Hillenbrand group I sou rces. The puffed-up inner 
rim in the Dullemond et al. (2001) model naturally explains 
this excess. In other words, no dynamically active model is 
needed to explain this bump. Furthermore, the outer parts of the 
disk can be flared, hence increasing the angle a under which 
the stellar light impinges onto the disk surface. The UV flux 
can penetrate deeper into the disk, thus increasing the geomet- 
rical thickness of the [O i] emission layer. As we have shown 
in the present paper, the observed [O i] 6300A emission profile 
can be explained as being due to the photodissociation of OH 
molecules and the subsequent non-thermal excitation of oxy- 
gen atoms in the atmosphere of a rotating flared disk. 



Appendix 

The fractional level population of two levels a and b of an atom 
in thermal equilibrium at a temperature T[K] is given by 



— = 8b 1 ex (-—) 
N a g fl i + ^ eXP l Tl 



(11) 



In this formula, g a and gb are the statistical weights, Aba is the 
Einstein transition rate for the transition between level b and a 
and Eba is the energy of the transition in K. The factor n is the 
number density of the collisional partner, y/ M is the collisional 
rate for the transition and the collisional partner. When more 
than one collisional partner is involved, the quantity nyt a is 
equal to the sum of these values for each interacting species i: 
ny ha = 2,n,rL- 



5 Hillenbra nd et alj have defined their groups based on the near-IR 
spectral slope. There is no direct link between both classifications: 
Hillenbrand group I and II contain Meeus group I as well as group II 
sources. 



22 



Acke, van den Ancker & Dullemond: [O i] emission in HAEBE disks. 



To determine the fraction F of neutral oxygen in the l Di 
state in the five-level model assuming thermal excitation, one 
applies 



(12) 



These ratios can easily be computed using equation dl It . the 
atomic data and the fact that 



Na 

N,, 



iNb 



(13) 
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